Long-range spin transport in superconductors 
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We report on non-local transport in multiterminal superconductor-ferromagnet structures, which 
were fabricated by means of e-beam lithography and shadow evaporation techniques. In the presence 
of a significant Zeeman splitting of the quasiparticle states, we find signatures of spin transport over 
distances of several fim, exceeding other length scales such as the coherence length, the normal- 
state spin-diffusion length, and the charge-imbalance length. The relaxation length of the spin 
signal shows a nearly linear increase with magnetic field, hinting at a freeze-out of relaxation by the 
Zeeman splitting. We propose that the relaxation length is given by the recombination length of 
the quasiparticles rather than a renormalized spin-diffusion length. 



The creation and and manipulation of spin-polarized 
currents forms the basis of spintronics applications [if. 
One key ingredient is the ability to transport spin cur- 
rents over mesoscopic length scales, which are usually 
limited by spin-flip or spin-orbit scattering processes. Su- 
perconductors are particularly interesting for spin injec- 
tion experiments due to the possibility to create almost 
100 % spin polarization enhanced spin relaxation 
times [3j], and the separation of spin and charge degrees 
of freedom Only few experiments on spin injection 
into superconductors have been reported so far (see [5|~0] 
and references therein). Both anomalously short 0] and 
anomalously long Q relaxation times as compared to the 
normal state have been reported. Here, we report on in- 
vestigations of spin transport in superconductors in the 
regime of large Zeeman splitting. In this regime, a cur- 
rent spin polarization of 100 % can be achieved 0, which 
is implicitely assumed in the classic experiments on spin 
polarized tunneling in high magnetic fields |8j . We study 
in detail the diffusion of spin-polarized quasiparticles by 
using a nonlocal detection with ferromagnetic electrodes 
as a function of contact distance, temperature and mag- 
netic field, and present evidence for spin transport over 
surprisingly long distances. 

Our samples were fabricated by e-beam lithography 
and shadow evaporation techniques. They consist of 
a thin superconducting (S) aluminum strip of thick- 
ness £ai ~ 10 — 15 nm, which was oxidized in situ to 
form an insulating (I) tunnel barrier before being over- 
lapped by several ferromagnetic (F) iron contacts (^Fe ~ 
15 — 25 nm). In addition a copper layer (icu ~ 30 nm) 
was evaporated under a third angle to reduce the resis- 
tance of the iron leads. Consistent results were obtained 
from nine samples of slightly different designs. We focus 
here on one sample (labeled FISIF) for which the most 
complete data set was recorded, and data from a ref- 
erence sample with normal-metal (N) copper contacts, 
but otherwise similar parameters (labeled NISIN). Fig- 
ure [T] shows a scanning electron microscopy images of the 




FIG. 1. (color online) Scanning electron microscopy image 
of a sample together with the measurement scheme with the 
injection (inj) and detection (det) circuits. 

FISIF sample, together with the experimental scheme. 
The sample has five contacts, spanning contact distances 
d from 0.5 to 8 /xm. 

The local (g\ oc = dl- ln j / dVi n j) and non-local (g n i = 
dld e t / dVinj) differential conductance for different contact 
pairs was measured by standard lock-in techniques in a 
setup described elsewhere 0, [l(|. Measurements were 
performed in the superconducting state of aluminum at 
temperatures down to T — 50 mK, and with an in-plane 
magnetic field B applied parallel to the ferromagnetic 
wires. For all data shown here the magnetization of the 
iron wires is aligned parallel to the magnetic field. We 
also performed non-local spin-valve experiments in the 
normal state at T — 4.2 K (not shown), from which the 
normal-state spin diffusion length An = 370 ± 10 nm 
and the spin polarization of the tunnel conductance 
P={G i - G t )/(G i + G f ) = 0.19 ± 0.05 were obtained. 
Here, Gf 4. are the junction conductances for each spin. 

Figure EJJa) shows the local differential conductance of 
one contact as a function of the injection bias voltage V5 n j 
for different applied magnetic fields B at T = 50 mK. For 
small B pronounced gap features at V « ±205 juV are ob- 
served as well as a negligible subgap conductance. Upon 
increasing the magnetic field, the gap features broaden 
due to orbital pair breaking, and for B > 0.5 T the Zee- 
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FIG. 2. (color online) (a) Local differential conductance 
ffioc = d/inj/dVinj of one junction as a function of injector bias 
Vinj for different applied magnetic fields B. (b) The same data 
plotted on a color scale. The lines indicate the regions were 
a single spin band dominates conductance. 



man splitting is seen. To describe our data, we have 
used the standard model of high-field tunneling Q to ob- 
tain the spin-dependent density of states n a (E), where 
a = ±1 stands for spin up and down, respectively. From 
n a (E), we calculate the current for each spin 
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where Gn = G^ + G^ is the junction conductance in the 
normal state, and fo is the Fermi function. The total 
charge current through the interface is then I = /-j- + Jj., 
and the spin current is proportional to I s = Jj- — /(.. 
Fits of this model to the measured conductance spectra 
yield Gn, the pair-breaking parameter T, and the spin- 
orbit scattering strength b so . Details of the fit procedure 
have been given previously [t| [13] • The spin polarization 
P obtained from these fits is consistent with the figure 
obtained from the spin-valve experiments. Figure [2jb) 
shows a contour plot of the complete dataset of the local 
conductance as a function of bias and magnetic field. In 
the wedge-shaped regions indicated by the lines, a single 
spin band dominates conductance. 

Next, we want to focus on the nonlocal differential con- 
ductance. To eliminate the effect of small variations of 
the junction conductances, we plot the normalized non- 
local conductance g n i/Gi n jGdet throughout the paper. In 
Fig- Eta) <7„i/Gi n jGdet is displayed as a function of the ap- 
plied bias voltage Vinj for different magnetic fields B and 
a contact distance d w 1 /jm. The data were measured 
simultaneously with the local conductance of Fig. EJa), 
in the configuration shown in Fig. [1] For comparison, 
we show data obtained from the NISIN reference sample 
in Fig. Elb). At B — 0, there is no conductance below 
the gap, and above the gap, both the FISIF and NISIN 
samples show a nearly linear increase due to charge im- 
balance [i]. With increasing magnetic field, the charge 
imbalance signal decreases, as clearly seen for the NISIN 
sample. The FISIF sample shows a qualitatively differ- 
ent behavior: (i) in the bias range corresponding to the 



FIG. 3. (color online) Normalized nonlocal differential con- 
as a function of injector bias VS n j for 
different applied magnetic fields B for one pair of contacts 
(a), nonlocal conductance of a pair of contacts of the NISIN 
reference sample (b), the data from panel (a) plotted on a 
color scale (c), and calculated differential spin current (d). 



Zeeman splitting, a positive peak arises for V m j < 0, 
and a negative peak for V; n j > 0; (ii) for higher bias 
I Vinj I > 300 /xV, an additional asymmetry evolves on top 
of the charge imbalance signal. While the observation 
(i) is systematic for all nine samples, (ii) was observed 
only in a few samples, whereas other samples showed the 
symmetric charge imbalance signal also seen in the NISIN 
sample at high bias. In the following we therefore only 
concentrate on the asymmetric peak features. Upon in- 
creasing the field, the peak heights increase gradually to 
their extremal values at B ~ 0.5 — 0.75 T, before the 
peaks start to decline, broaden and move inwards, simul- 
taneously. The positive peak (at negative bias) is slightly 
larger than the negative peak (at positive bias). Above 
the critical field B c « 2.15 T the asymmetric features 
disappear and one finds a small bias-independent signal 
(not shown). 

Figure [UJc) represents the same data as in Fig.lUJa) as 
a contour plot with additional magnetic field steps. The 
lines are duplicated from Fig. [U[b). As can be seen, the 
asymmetric conductance features are limited almost en- 
tirely to the magnetic-field and bias region of the Zeeman 
splitting of the quasiparticle density of states. For com- 
parison, in Fig. EJd) we show the differential spin current 
dI s /dVi n j calculated using ([T]) and parameters obtained 
from the fits to the local conductance of the injector. 
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FIG. 4. (color online) (a) Normalized nonlocal differential 
conductance <J n i/Gi n jGdet as a function of the bias voltage 
VSnj for different contact distances d at fixed magnetic field B 
(b) peak area A (for the peak at Vi n j < 0) as a function of B 
for different d (c) semi- logarithmic plot of A as a function of 
d for different B, together with fits to an exponential decay 
(d) relaxation length As derived from these fits as a function 
of B, together with the charge imbalance length Aq* of the 
reference NISIN structure, and the normal-state spin diffusion 
length An derived from non-local spin-valve experiments. 



Fig. 2][a) shows the normalized nonlocal conductance 
as a function of the bias voltage for different contact 
distances d and one fixed magnetic field and tempera- 
ture (B = 1 T and T = 50 mK). In general, increas- 
ing the distance between injector and detector only de- 
creases the amplitude of the signal, but does not affect its 
overall shape. Only the signal for the shortest distance 
d = 0.5 /im slightly deviates in shape. 

To further characterize the observed peaks, we have 
analyzed the heights and areas of both peaks as a func- 
tion of magnetic field, temperature and contact distance. 
Height and area for both peaks yielded similar results, 
and we only show the peak area A of the positive peak 
in the following. Fig. 0Jb) shows A as a function of the 
applied field B for different distances d at T — 50 mK. 
A rises monotonically with increasing magnetic field, un- 
til at B r; 1.5 — 1.75 T it reaches a maximum and then 
rapidly goes back to zero. Again, this behavior is similar 
for all distances, with slight deviations at d = 0.5 /xm. 

Fig. |4jc) shows the peak area as a function of contact 
distance on a semi-logarithmic scale for different mag- 
netic fields. The solid lines are fits to an exponential 
decay, from which we obtain the relaxation length As- 
As can be seen, the quality of the fits is very good up to 
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TABLE I. Four terms in the tunnel Hamiltonian involving the 
state kcr in the superconductor. 



B rs 1.5 T, except for the shortest distance d = 0.5 /xm, 
which was therefore excluded from the fits. At fields 
above 1.5 T, the quality of the fits declined (not shown). 

Fig. IDJd) shows the relaxation length As as a function 
of the applied magnetic field B for two different tem- 
peratures. It is found to be about 2 [im for small mag- 
netic fields and then increases nearly linearly with the 
applied magnetic field up to around 7 /xm for B s=s 1.5 T. 
At higher fields, As appears to decline, but as men- 
tioned above the fits were not very good in this region. 
As can be seen, As is almost independent of tempera- 
ture up to 350 mK. For comparison, we also show the 
charge-imbalance relaxation length Aq» obtained from 
the NISIN reference sample, and the spin-diffusion length 
An from the normal-state spin-valve experiments. 

For a simple qualitative model, we consider a BCS su- 
perconductor in high magnetic fields, including the effect 



of the Zeeman splitting 11], but neglecting orbital pair- 
breaking and spin-orbit scattering for simplicity. In the 
normal state, the electron energy relative to the chem- 
ical potential \i is given by ek = h 2 k 2 /2m — /i. The 
quasiparticle energies in the superconducting state are 
Ek<r = Ek + a-fXsB, where = \Je^ + A 2 . Tunneling is 
described by a straightforward extension of 0, HH to the 
case of finite Zeeman splitting: The tunnel Hamiltonian 
for spin-conserving tunneling between states kcr and qrr 
in the superconductor and ferromagnet, respectively, is 



H T = ^2 T qk(T c k(T c q(T + H.c 

qk(j 



(2) 



Electron operators c in the superconductor are written 
in terms of the quasiparticle and Cooper pair operators 
7 and S as 



(3) 



where k=— k, a — ~ a, and the coherence factors are 



« = 1 



The four terms appearing in the tunnel Hamiltonian for 
the state kcr in the superconductor are listed in table U 
where /ko- is the quasiparticle distribution function in the 
superconductor, and /o is the Fermi function describing 
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occupation in the ferromagnet. Summing up all contri- 
butions to the current for both spin directions yields 

1 = - E / (/- " MEea ~ eV)) 

(7 ( 5 ) 

-G w v 2 AU<y~ k{Eea+eV))} de, 

where we have replaced sums over k by integration over 
e. We now group terms with equal energy E inside 
(e < 0) and outside (e > 0) the Fermi surface, and intro- 
duce transverse and longitudinal distribution functions 
/iJ ) = (U ~ fea)/2 and f£ } = (/„ + jW))/2 in the su- 
perconductor. Setting eV = for the detector junction, 
we obtain 

1 = ~ Y. + G °) Qa + ( G ° ~ G °) S °\ > ( 6 ) 

cr 

where 

Q*=/ («S-«?)/^«fc (7) 

Je>0 

S*= [ - fo(Ee*j) de (8) 

Je>0 y J 

describe the charge imbalance and spin accumulation in 
each band. The first term in ([6]) is the usual charge- 
imbalance signal observed in both the NISIN and FISIF 
sample. The second term is non-zero only for a ferro- 
magnetic detector junction and gives rise to the observed 
peaks, as explained below. 

Current injection in the bias regime of the Zeeman 
splitting creates both a charge imbalance QT and spin 
accumulation in the lower Zeeman band, while the 
upper Zeeman band remains unoccupied. After diffusion 
over a distance d > Xq* along the superconductor we 
have Qt ~ 0, but there is still a finite spin accumulation 
Si > (spin-charge separation). The detector current is 
then 

/det«^(G t -G^)% (9) 

We would now like to discuss the salient features of our 
experiment and some open questions within this model, 
(i) We expect S± oc |ii n j|j since for both bias directions 
there is an excess population of quasiparticles. Conse- 
quently, the sign of the detector current depends only on 
the sign of the spin polarization P. Since both Sj, and 
Jdet are even functions of bias, the differential signals dis- 
played in Figs. [3] and |4] are odd functions, explaining the 
asymmetric peaks, (ii) Once the upper Zeeman band is 
reached by the bias voltage, it yields an opposite con- 
tribution, canceling the signal. Therefore, the signal is 
restricted to the region of the Zeeman splitting, as seen 
in Fig. [3jc). (hi) The peak height is proportional to the 
injector current. Therefore, the peak at negative bias is 
slightly larger than the peak at positive bias. The peak 



height ratio should be (1+P)/(1— P) « 1.5, which is con- 
sistent with the data shown in Figs.[3]a) and|4|a). (iv) In 
the bias window of the Zeeman splitting, the relaxation of 
the quasiparticles is strongly suppressed. They can either 
diffuse out of the wire into the superconducting reser- 
voirs, or recombine to Cooper pairs. This is qualitatively 
different from a renormalization of the normal-state spin- 
diffusion length due to coherence factors [3]. For recom- 
bination, the quasiparticles need a spin-reversed partner 
from the upper Zeeman band, which is almost unoccu- 
pied. Consequently, As is very large. Weak relaxation 
paths might come from spin-orbit scattering [l3[ or from 
inelastic spin flips to the upper Zeeman band. The lat- 
ter would also explain the increase of As with magnetic 
field, since larger energy transfer would be needed for in- 
creasing Zeeman splitting. A quantitative model for the 
relaxation mechanism is an open question to theory, (v) 
Asymmetric peaks in the nonlocal conductance were also 
predicted for the competition of crossed Andreev reflec- 
tion and elastic cotunncling in FISIF structures in the 
presence of Andreev bound stat es fl3 . [l5| . Bound states 
generated at the FIS interfaces [lOj], however, should ex- 
tend no further than the coherence length into the super- 
conductor, which is clearly inconsistent with the length 
scales observed in our experiment. Nevertheless, a com- 
mon feature of bound states and Zeeman splitting is the 
existence of a spin-polarized density of states in the su- 
perconductor, and the conductance features look intrigu- 
ingly similar. Another question to theory is therefore 
whether a coherent picture of both effects can be ob- 
tained. 

In conclusion, we have demonstrated spin injection and 
transport in superconductors in the regime of large Zee- 
man splitting of the density of states. We have found 
an asymmetric nonlocal conductance signal that can be 
modeled by spin accumulation in the lower Zeeman band 
of the quasiparticle dispersion. The relaxation length of 
the spin signal exceeds the normal-state spin-diffusion 
length by one order of magnetitude, and we propose that 
it is determined by the recombination rate of the quasi- 
particles. 
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